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I discuss the quark mass dependence of various baryon properties derived from chiral perturbation 
theory. Such representations can eventually be used as chiral extrapolation functions when lattice 
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in loop and contact term contributions, the latter given in terms of low-energy constants. I stress 
the importance of utilizing phenomenological input to constrain a certain class of low-energy 
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splitting is counted as a small parameter. I also discuss briefly the pion mass dependence of the 
nuclear force as derived from chiral nuclear effective field theory. 
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1. Approaching the chiral limit in QCD 

Chiral perturbation theory (CHPT) is the appropriate tool to analyze the quark mass depen- 
dence of baryon properties. CHPT is the effective field theory (EFT) of QCD, based on its sym- 
metries and their realization [jj], ^, f3|. A central role in QCD is played by the spontaneously and 
explicitly broken chiral symmetry for the three light quark flavors up, down and strange. The 
vacuum is not invariant under this global SU(3)lxSU(3)# symmetry and the spectrum therefore 
contains eight almost massless pseudoscalar Goldstone modes. These are not exactly massless 
because the quark masses are non-zero but have very small values compared to typical strong in- 
teraction scales. The consequences of this explicit and spontaneous chiral symmetry breaking can 
be analyzed in the framework of CHPT. S -matrix elements and transition currents are expanded in 
small parameters - external momenta and quark (meson) masses with respect to the chiral symme- 
try breaking scale A^ ~ 1 GeV. Matter fields (that is all particles that are not Goldstone bosons) 
can also be considered in CHPT. A special role is played by the ground state baryons - they can 
not decay strongly and thus one is able to consistently deal with the new mass scale related to these 
fields - the chiral limit value of the baryon mass is neither non-vanishing nor small compared to 
the scale of chiral symmetry breaking A%. Resonances, may they be meson or baryons, can not 
be incorporated so directly, one either has to resort to resummation schemes or consider additional 
small parameters that do not vanish for vanishing quark masses in the construction of the effective 
field theory. For reviews on CHPT including also matter fields I refer to [Q, |5[ |6[ 0, H and the 
foundations of CHPT are discussed in depth in [g]. 

We now consider the approach to the chiral limit. For the following discussion it is important 
to differentiate between the SU(2) and the SU(3) chiral limit (c.l.) of QCD. More precisely, in the 
SU(2) c.l. we have m u ^md — > with m s fixed at its physical value and the heavy quarks decoupled; 
while in the SU(3) c.l. m u ,md,m s — ► with the heavy quarks decoupled. There are a few general 
observations about the chiral limit of QCD. First, S-matrix elements exist in the chiral limit for 
arbitrary momenta. Second, as it is well-established, the approach to the chiral limit is non-analytic 
in the quark masses m q , leading to the famous chiral logarithms and other non-analytic behavior 
of certain observables. Third, in the fundamental paper [ |10| ] a decoupling theorem was derived. 
It states that the leading chiral non-analytic terms stem from pion (Goldstone boson) one-loop 
graphs coupled to pions (Goldstone bosons) or nucleons (ground state baryons). This theorem has 
various immediate consequences: resonances like the p, the A, . . . decouple to leading order. This 
puts severe constraints on the construction of EFTs with resonance fields. Related to this is the 
observation that chiral limit of QCD and the large N c limit do not commute - we will come back to 
these issues later. 

2. One-loop representation of baryon observables 

We now consider a complete one-loop (fourth order) SU(2) calculation. Any nucleon observ- 
able - given in terms of tree and loop diagrams - takes the form 

ff(M n ,F n ;mo,g A ;ci,di,ei) , (2.1) 

where some of the quantities entering this expression are themselves quark mass dependent. Specif- 
ically, the pion mass and the pion decay constant are given by M\ = (m„ + m^Bo + @{m 2 q ) and 
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F% = Fo [1 + G(m.q)\, where Bo is related to the quark condensate and Fq is the pion decay constant 
in the SU(2) chiral limit. Further, mo is the nucleon mass in the SU(2) c.l. and gA is the leading 
order (dimension one) axial-vector pion-nucleon coupling, gA = go [1 + 0(m q )]. Also, the Ci,di,et 
are low-energy constants (LECs) of the dimension two, three, and four chiral effective pion-nucleon 
Lagrangian if eff , 

■^eff = ~^%N ~^%N -^TtN + • • ■ ) (2-2) 

where the superscript (n) denotes the chiral dimension (the purely pionic part of the Lagrangian is 



not made explicit here). The complete effective Lagrangian is given in fllip. For later use, I display 
the dimension one and two terms: 



=w(i0-mo + ±goi/Y 5 )y 



7 



1 , v _„„ . s 1 



= £*c,-Op ) V = * (c 1 (x + )+c 2 {--^(u ll u v )D^ + h.c.) +c 3 Uu-u) 
+ c A l - [b m , uy] + c 5 x+ + c 6 -^F^ V + a <^ + v> ^ V V , 

23 118 

= E**0f , JSf$ = £^C> (4) ^ , (2.3) 

i=\ 1=1 

with f/(x) = w 2 (x) collects the pion fields %{x) = x a ■ % a {x), *P(jc) the isodoublet nucleon field, 
U[i = iu^DpJJifi = —id^n/Fo + ... and the chiral covariant derivative. The Of are monomials 
of chiral dimension n. Furthermore, x+ parameterizes the external scalar source including the 
quark masses, %+ ~ 2Z?odiag {m u ,m£) + F+ v is the chiral and gauge covariant photon field 
strength tensor, and ( ) denotes the trace in flavor space (for more details, see [g 0). The LECs 
can be grouped into two distinct classes. Class I are the so-called dynamical LECs ~ <??,<9„,... 
which are non-vanishing in the chiral limit. Many of these are accessible via phenomenology, 
that is through the analysis of scattering data or decay processes. Class II LECs, the so-called 
symmetry breakers ~ m q ,ml,m c ,d^, . . ., vanish in the c.l. and they can best be obtained by varying 



the quark masses, i.e. in lattice simulations (see also the discussion in [12]). Note that the mixed 
LECs parameterizing operators with quark mass insertions and derivatives are counted as symmetry 
breakers because at fixed pion mass these can be absorbed in the values of certain dynamical LECs 
(with the exception of the leading mixed LECs at dimension three). The dimension two effective 



Lagrangian in Eq. ( |2.3| ) has two symmetry breakers and five dynamical LECs, from the latter two 
can only be determined in the presence of external fields. At dimension three, we have 4 symmetry 
breakers and 19 dynamical LECs - in fact, the symmetry breakers are necessarily of the mixed 
type and thus are easier available from phenomenological analysis as most other LECs of that type. 
Finally, at dimension four, we have 25 symmetry breakers (most of them of mixed type) and 93 



dynamical LECs (most of them containing external fields). For more details, see [11]. As we will 
see in the following, while the total number of LECs increases rapidly, this is not the case for most 
observables - this will become evident when I discuss the nucleon mass, the isovector magnetic 
moment or the axial-vector coupling constant. 
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3. Analysis of baryon properties 

Chiral perturbation theory supplies the quark mass expansion of any given observable to the 
required accuracy - provided the expansion parameter is sufficiently small. There is an impor- 
tant interplay between chiral loops and the LECs parameterizing the contact interactions. Mostly, 
one employs dimensional regularization to perform the necessary renormalization of the loop con- 
tributions in a symmetry-preserving manner (although other schemes are viable, too). While an 
observable & does not depend on the scale of dimensional regularization A, that is d(7 '/dX = 0, 
this is in general not the case for the individual contributions from the loops or the counter terms. 



An example of this is the isovector charge radius of the nucleon as discussed in Q13[]. It is therefore 
absolutely necessary to include all terms in harmony with the underlying symmetries, and not just 
the chiral logs or other non-analytic terms. Also, it has become evident over the years that in baryon 
CHPT one has to perform complete one-loop (fourth order) calculations to obtain a precise enough 
representation at physical pion/kaon masses (in SU(3), this statement is not true in general). For 
reviews, see p4| , [15[ |. 

Before continuing, let me make the following disclaimer . I will discuss continuum extrapola- 
tions, but no finite lattice spacing or finite volume effects Also, I will not discuss (partial) quench- 
ing, that means all lattice data are taken as full QCD. Further, I will be mostly interested in the 
range of applicability (theoretical uncertainty) rather than precision fits - for reasons that become 
apparent in the following - and finally, I eschew models here - end of disclaimer . 

3.1 Quark mass expansion of the nucleon mass in SU(2) and SU(3) 

Consider the expansion of the nucleon mass in the light quark mass m = m u = m c i with m s 
fixed. By means of two-flavor CHPT, this is mapped onto an expansion in the pion mass. Note that 
we consider strong isospin violation (m u ^ m c {) later. The fourth order calculation gives [ |I6| ] 

m N = m Q - 4 Cl Ml - -f±-§ + * 1 M> / -^ + k 2 M 4 n + 0{M 5 % ) , (3. 1) 

with mo the nucleon mass in the chiral SU(2) limit and ki and k 2 are combinations of second 
(ci,C2,C3) and fourth order LECs (<?i), k\ = — (3/32n 2 F 2 )(— 8ci + c 2 +4c3 +g\/mfj) and k 2 = 
—Ae\ + (3/\28n 2 F 2 )(c 2 — 1^ A jm^), respectively. The scale of dimensional regularization was set 



equal to the nucleon mass, A = m^. For simplicity, I have neglected in Eq.(3.1) the difference 



between the physical value of the pion mass and its leading term in the quark mass expansion, 



M 2 = M 2 [l + 0(m q )] (for details, see e.g. [17]). Note also that the dimension two LECs are 



scale-independent in any mass-independent regularization scheme, while the fourth order LEC is 
scale-dependent, e\ = e\{X). It is important to realize that the LECs ci,c 2 ,c?,, and e\ are strongly 
constrained from pion-nucleon scattering and to some extent form peripheral nucleon-nucleon 
scattering. More specifically, the chiral expansion of the %N amplitude converges best inside the 
Mandelstam triangle, however, to determine the values of the LECs, one must compare to a dis- 
persive analysis of the scattering data. This was done in []18[]. Alternatively, one can use the chiral 
representation of the amplitude in the physical region close to threshold, here one is faced with 
certain inconsistencies in the data basis. The most detailed fits based on these data were performed 
in [p^]. One can also pick out certain threshold parameters that are particularly sensitive to certain 
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LECs, as it was done in Q2(J|. In addition, the peripheral partial waves in elastic nucleon-nucleon 
are sensitive to the two-pion exchange and thus to some of the c,-. In this case, one has a large data 

- which, 



basis but must subtract the dominant one-pion exchange, see [ |21| , [22[ ] (see also Ref. [|2 
however, does not enter my average). Putting all this information together, we obtain for the c, 
(i = 1, . . . ,4) (all values in GeV" 1 ) 



-0.9 



t-0.5 
-0.2 > 



C2 



3.3 ±0.2 , c 3 = -4.7 



+1.2 
-1.0 ' 



C4 = 3.5 



+0.5 
-0.2 > 



■2ci + — + c 3 



-2.1 



+ 1.3 



L -l.l • 
(3.2) 

A few remarks are in order. First, the symmetry breaker c\ and the isoscalar combination (last 
number) are least accurately determined - as expected. Also, the LEC e\ is only badly determined 
from 7lN — > 7iN in the threshold region, the analysis of [|4j] gives a negative value of natural size. 
It is also important to notice that the numerical values of the LECs c, can be well understood in 
terms of resonance saturation (baryon and meson resonance excitations), in particular, the fairly 
large values of C2,3 (£4) are generated mostly by A(1232) (p) exchange [ffij). Before discussing the 



quark mass dependence of the nucleon mass as given by Eq. ( p.lD , it is important to point out that 
for the physical pion mass, the higher order corrections are small. More precisely, for the central 
values of the LECs given in Eq. (3_^2) (and using e\ = — 1 GeV~ 3 ), the contribution quadratic, cubic 
and quartic in M n amounts to 70, —17 and — 4MeV, respectively. It was even shown in [25] that 



the fifth order corrections are tiny (at the physical point). I will return to this issue of convergence 
later. 



a) 
O 




0.4 

M [GeV] 



0.2 0.3 0.4 

M [GeV] 



Figure 1: Pion mass dependence of the nucleon mass. Left panel: SU(2) analysis. The blue solid line refers 
to the best fit and the dashed black lines to the theoretical uncertainty as discussed in the text. Red circles: 
MILC data [Q, black squares: CP-PACS data [Q. Right panel: SU(3) analysis. The red line gives the best 
fit as discussed in the text and the black lines give the theoretical uncertainty under the constraint that the 
nucleon mass takes its physical value for the physical quark masses. Green squares: MILC 2001 data 
purple triangles: MILC 2004 data 



In the left panel of Fig. | we show the results from Ref. [£3]]. A best fit to the MILC 2001 @ 
and the CP-PACS data [[T7]] is obtained with a = -0.9, c 2 = 3.2, c 3 = -3.5 (all in GeV" 1 ) and 
e\{m^) = — lGeV~ 3 . Note that this is not a chiral extrapolation in the true sense of the word 
because the constraint that the physical value of is obtained at the physical pion mass was 
imposed. Note also that the MILC data are from an SU(3) simulation, however, they have also 
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given an SU(2) reference value for one pion mass that differs only little from the corresponding 
SU(3) result. The resulting SU(2) c.l. value of the nucleon mass is mo ~ 0.88 GeV. The dashed 
lines in the figure reflect the theoretical uncertainty from the variation of the LECs, consequently, 
we obtain a moderate/large theoretical uncertainty for M n above 400/550 MeV. Clearly, at higher 
pion masses the higher order terms become too important for such fits to make sense. For related 



work on the nucleon mass in SU(2), see [29, BQ]. The theoretical uncertainty in the determination 



of the pion mass dependence of the nucleon mass is also discussed in p 1Q- 

Next, I consider the extension to the three-flavor case. Before discussing the nucleon mass, 
some general remarks are in order. As compared to the SU(2) analysis, one has more fields and 
operator structures in the EFT, and consequently, more LECs. As examples, consider the leading 
order axial meson-baryon coupling in SU(2) and SU(3) 

{gAW^f¥^D(BY^f{u^B})+F(BY^Y 5 [u^B}) , (3.3) 

or the leading symmetry breakers 

ci\jr(x+W->b (BB)+b D (B{x + ,B})+b F {B\z + ,B]) , (3.4) 

with \P the nucleon doublet and B the baryon octet coupled to the pion triplet and the Goldstone 
boson octet, respectively. Note that the LEC bo can be absorbed in the value of the octet baryon 
mass in the chiral limit, but need to be kept separately when one also analyses the pion- and kaon- 
nucleon sigma terms [p2|]. These various operators in SU(2) and SU(3) are related by matching 
conditions - in the SU(2) case the strange quark effects are buried in certain local operators. Such 
matching relations have been worked out in [p3jj, to leading order, one finds, e.g. 



nucleon mass in the SU(3) chiral limit : mo = mo [1 + &(m s )] , 
leading axial coupling : gA = D+F+ &(m s ) , (3.5) 

leading symmetry breaker : c\ = bo + j(bD + bp) + 



/mo 



2 

The full matching conditions to fourth order in the chiral expansion are given in p3[]. These are 
important constraints that should be implemented in any SU(3) analysis. I return to the nucleon 
mass. Ref. [ |33| ] contains the complete fourth order expressions (including isospin breaking terms) 
for the baryon octet masses and sigma terms. Based on that work, in [|4|] fits to the MILC 2001 
data based on two different regularization schemes (cut-off and dimensional regularization) were 
performed. Again, the constraint to obtain the physical value of the nucleon mass at the physical 
quark masses was imposed. The dimension two LECs called bj were taken from the earlier work 



in [35] and the three combinations of the dimension four LECs (called dj) were varied to obtain a 



best description of the three-flavor MILC 2001 data, see the red line in the right panel of Fig. [Tj In 



the figure are also shown the MILC 2004 data [ 28 ] at lower quark masses - these can not be fitted 



with the constraint (as already noted by the MILC collaboration using a simplified mass formula). 
Taking into account also the kaon mass dependence of the nucleon mass and the uncertainty due to 
the MILC 2004 data, we obtain the following ranges for various (isoscalar) quantities: 

SU(3) c.l. value of m N : 710 MeV < m < 1070 MeV , 

Pion - nucleon sigma term : 39.5 MeV < 0^(0) < 46.7 MeV , (3.6) 

Strangeness fraction : 0.07 < y < 0.22 , 
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physical value 
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[GeV] 



0.6 
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■ MILC01 

▲ MILC04 



0.2 



0.4 

M [GeV] 



0.6 



Figure 2: Pion mass dependence of the A mass (left panel) and of the E mass (right panel). Green 
boxes/purple triangles: MILC 2001/2004 data. The blue circle denotes the physical value of the corre- 
sponding baryon mass. For further details, see the text. 



with Gjin(0) = (N(p)\m(uu + dd)\N(p)} (with N(p) a nucleon state of momentum p) and the 
strangeness fraction in the proton, y = 2{p\ss\p) / (p\uu + dd\p) , can be obtained from (7^(0) = 



°b/(l — y) an d using Oq = (37 ± 6)MeV from H33]. It was concluded in [34] that for pion masses 



less than 400 MeV the convergence of the chiral expansion is fine and that it is acceptable for 



masses below 550 MeV. This agrees with the findings of the SU(2) analysis in [13]. 



3.2 Quark mass expansion of the baryon octet masses 

We now consider the other octet members with strangeness (A, £, E) based on the work of Ref . [Q] . 
In that paper, the fourth order LECs were fitted to the nucleon mass only, that means the Gell- 
Mann-Okubo relation for the baryon masses was not imposed. Consequently, some of the masses 
come out somewhat off their empirical values. More precisely, while the £ mass is well reproduced, 
the A and E masses come out by about 10 — 15% too high. To get a handle on the theoretical 



accuracy, we also use the values for the d\ from Q35[], in that case all masses are exactly reproduced. 



The resulting prediction for the pion mass dependence of the A and the S in comparison to the 
MILC data are shown in Fig. |[ The red/black lines refer to the optimal set of the LECs from the 
nucleon mass fit/to the LECs from [[3^]. We note in particular that the pion mass dependence for 
the S is much flatter as one would expect from the MILC data. This is not unexpected - the S only 
contains one valence light quark and should thus be less sensitive to variations in the pion mass. 
Clearly, one could improve this description by fitting directly to these particles. Still, the pion mass 
dependence of the S as given by the MILC data is a mystery to be resolved. 

So far, I have considered the isospin limit m u = mj. I briefly discuss the implications of 
unequal light quark masses for the strong neutron-proton mass difference. In the SU(2) case, it is 
given to leading order by 

(m n - m p f mn z = 4c 5 B (m u - m d ) + 0((m u - m d f) (3.7) 

with C5 the dimension two LEC accompanying the operator %+ = %+ — h{X+) ~ #o( m « — nid), see 



also Eq. (2.3). Utilizing the Cottingham sum rule to estimate the electromagnetic splitting, one 
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finds pt] 



-0.09 ± 0.01 GeV 



(3.8) 



The small value of this LEC is a reflection of the suppression of isospin violation in QCD, which 
scales as (m u — m ( i)/A x . (Note that Weinberg obtained a slightly larger value in his seminal paper 



because he employed SU(3) relations [Bq]). The fourth order corrections to Eq. (3.7) have been 



worked out, see e.g. Q3Vp . In SU(3), the expression for the strong neutron-proton mass is much 
more complicated, since one has to treat 71° — r\ and A — Z° mixing in terms of the parameter 
tan2e = (v3/2) (m u — m^) / (m s — m). Again, one can derive a matching condition [p3|], to leading 
order it reads (for related work, see the pioneering analysis in [[38]] and also [39]) 



c 5 = bp + bp- + ff{m S) m s \nm s ) 



(3.9) 



Needless to say that for obtaining precise mass splittings, one needs to consider also the electro- 
magnetic corrections (for an early analysis in the meson sector, see [fl(!|]). 

3.3 Quark mass expansion of other nucleon properties 

Besides the nucleon mass and sigma term, chiral extrapolations for various other nucleon 
observables based on chiral perturbation theory (or extensions thereof) have been considered in 



the literature, e.g. nucleon magnetic moments and the electromagnetic radii [ |41| , |42| , |43| ] or the 



iso vector axial-coupling constant ||44j, |45|, |4g]. I will not review these in detail but rather make 



some more general remarks on the present status of extrapolation functions derived from baryon 
chiral perturbation theory. There are two important issues which require special attention: 
|~T| Given the scarcity of data at low pion masses (say below M K ~ 400 MeV), one should per- 
form global fits to a variety of observables at sufficiently small quark masses. The important 
point is that the LECs - the parameters of these fits - relate many observables, they are the 
same for all processes and can therefore not be determined independently. Furthermore, one 
should incorporate as much phenomenological input as possible, in particular for the dynam- 
ical LECs, see the discussion of the nucleon mass in the preceeding sections. Needless to say 
that one must be in a regime where higher order terms stay sufficiently small. Attempts to 



® 



go beyond this regime are necessarily model-dependent, see e.g. Q47Q and references therein. 
Results should be independent of the regularization scheme, since these can differ only 
by higher order terms if symmetries are properly implemented. For example, the so-called 



infrared regularization (IR) of baryon CHPT [17] resums all kinetic energy insertions ~ 
p 2 /2mN whereas the latter are expanded order-by-order in the heavy baryon approach (there 
are a few circumstances for which the heavy baryon approach does not converge, but these 
are not of relevance here). Similarly, cut-off schemes or finite-range-regulators as employed 



by the Adelaide group, see e.g. [|48[], resum certain higher order corrections (if implemented 
properly). For small pion masses as defined above, all these schemes must give the same 
result. 

Let me illustrate these issues on the example of the nucleon isovector anomalous magnetic moment 
K v . It has been worked out to fourth order in [|^] (the pertinent Feynman diagrams are shown in 
the left panel of Fig. ||) 



K v = c 6 - l6m N M n em - 



..loop, 3 , „loop.4 



(3.10) 
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(a) 



(d) 



(g) 



(b) 



(e) 



(h) 



(c) 



(0 



(i) 



10 



8- 













physical value 


• 


QCDSF 




Figure 3: Left panel: Second (a), third (b) and fourth (c-i) order diagrams contributing to K v , The 
solid/dashed/wiggly lines stand for nucleons/pions/photons, respectively. The • (■) denotes a dimension 
two (four) insertion. The x denotes an dimension two insertion with fixed coefficients. Right panel: Pion 
mass dependence of K v as described in the text: The upper (lower) line corresponds to the c.l. (physical) 
values of the pion decay constant and the nucleon mass. The diamond represents the physical value of K v . 
The data are from QCDSF [ph. 



where I have made explicit the contributions from the dimension two ~ C6 (diagram (a) in Fig. ||) 
and four ~ contact interactions (diagram (c) in Fig. |3h. Here, C(, is nothing but the anoma- 
lous magnetic moment in the chiral limit, k£ = C(,. The fourth order loop contribution depends 
explicitely on the dimension two LECs ci_ (the tadpole graph (g) in Fig. |J), C4 (the tadpole graph 
(d) in Fig. [3b, and C(, (diagrams (e) and (f) in Fig. |3h. Two of these are constrained from pion- 
nucleon scattering, cf. Eq. fl3.2| ). The expression for K v contains the nucleon mass and the 
pion decay constant F n (in the loop contributions). These quantities depend themselves on the pion 
mass, niN = m^{M ]l ) and F n = F K (M 7[ ). To the order we are working, one can always replace the 
parameters in the effective Lagrangian (mo and Fo, the pertinent c.l. values, see e.g. Eq. ( |2.3[ )) by 
their physical values, the differences being of higher order in the chiral expansion. Thus, replacing 
{mo,Fo} by {m^^n} in the formula Eq. ( 3.10| ) is a way of investigating the convergence of the 
series, and I will do this here for fe v . The explicit expression for the pion mass expansion of 
is given in Eq. (3J_) and F K to GiM^ can be taken from [5C] (see also Ref. [ [51] ] and references 
therein) 



F 



l+X 



3L? 
4 



7L. 



4£ 2 

- ! - —h-^f+U 

3 iz/ z iz j ivz 

(3.11) 

where L = log(A 2 /^), X = MI/{\6k 2 F$) and r F (fx) is a combination of dimension six LECs. I 
use X = 0.5 GeV and r F (A = 0.5 GeV) = 3 in the following. The NLO LECs l t are tabulated in [£(]]. 
Guided by the trend of the QCDSF data [43] and imposing the constraint that K v = 3.71 [n.m.] at 
the physical pion mass, we set C(, = 5 and eio6 = 0.45. The resulting theoretical uncertainty due to 
the two choices for the pion decay constant and the nucleon mass are shown in Fig. 0. As before, 



-) 

12/ 



+ 



12 



13 



192 



+ r F (A) 



I am grateful to Veronique Bernard for providing me with these results. 
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below pion masses of about 400 MeV, the theoretical uncertainty is modest, but above it quickly 
increases - so that a model-independent extraction of K v using the particular data shown in the 
figure is not possible. 

Let me now consider the axial- vector coupling gA, which was much discussed in the recent lit- 
erature, see e.g. [44, 45, fjq ]. I present here some first results of an on-going study with Veronique 



Bernard Q52[] . There are two main issues concerning this particular observable. First, to my opin- 
ion, the existing data are at too high pion masses to draw definite conclusions about the range of 
applicability of the chiral extrapolation functions or the role of explicit spin-3/2 degrees of free- 
dom. Second, there is a much more direct problem with the pion mass dependence of gA . It has 



been worked out to fourth order first in [53]. Using the standard form of the effective pion-nucleon 



Lagrangian (see Eq.(2.3) and Ref. [11]) it reads 



gA = g0 



1 + - 



4M 



go 



igo + g , M n 



16tt 2 F ( 2 



In- 



16tt 2 F 2 



+ 



Mi 



24 K Fq mo 



(3 + 3gl - 4m c 3 + 8m c 4 ) 



(3.12) 



The dimension three coupling constant d^ can be determined from the process 7iN — > nnN [Q], 
however, not very accurately. For a typical value of the LEC die = — 1.76GeV~ 2 , the correction 

(2) 

quadratic in the pion mass is Ag A ~ 0.16, which is of typical size for an SU(2) correction (note 
that d\(, is the renormalized value of d\(, at X = M n , see fll9|]). It was, however, already pointed 
out in [53] that the corrections ~ M\ are already quite sizeable for the physical pion mass. Using 

(3) ■— -, 

the central values of the LECs C3 and C4, one finds Ag A ~ 0.32, using Fq = 86.5 MeV pOJ, mq = 
880 MeV [|l3|] and go = 1.2. This is a fairly large correction for an SU(2) quantity. It can be traced 
back mostly to the particularly large combination of the dimension two LECs —Amoc-i + %moC\ ~ 
41, and thus one can not expect this representation to be very accurate. This is further corroborated 
in Fig. |], where the pion mass dependence of gA is shown for various values of the LEC d\& and 
two values of go = 1.0 and go = 1.3, respectively. The strong increase of gA(Mx) with growing M n 
is due to the large correction cubic in the pion mass - this clearly narrows the window for applying 



Eq.(3.12) as a chiral extrapolation functions to even smaller pion masses than it is the case for the 
nucleon mass or the isovector magnetic moment discussed earlier. It should also be stressed that 
using an extrapolation function based only on the third order calculation (i.e. including the terms 
up to quadratic order in M n ) is insufficient - one does not even capture the largest correction at 
the physical value of the pion mass. On the other hand, it is conceivable that even higher order 
corrections are again of natural size - this, however, needs to be checked by an explicit calculation. 



A more detailed discussion of these issues will be given in [52] (see also [M4J) 



3.4 Quark mass expansion of the A mass 

The A(1232) is the most important baryon resonance. It is almost degenerate in mass with 
the nucleon and couples strongly to pions, nucleons and photons. It was therefore argued early 



that spin-3/2 (decuplet) states should be included in baryon chiral perturbation theory [55]. In that 
paper and subsequent works use was made of the heavy baryon approach, which treats the baryons 
as static sources. However, special care has to be taken about the decoupling of resonances in the 
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Figure 4: Pion mass dependence of the axial-vector coupling for various values of the LEC d\§ = 
—0.92, — 1.72,— 3.4GeV~ 2 and the chiral limit value of gA, go = 1.0 and 1.3, respectively. The filled or- 
ange circle denotes the physical value of gA at the physical pion mass. 



chiral limit, as discussed in Sec. [I] - this is frequently ignored in literature. This approach can be 
systematized by counting the nucleon-delta mass splitting as an additional small parameter, so one 
expands in the generic small parameter e []56|], with 

[ A x A x A x J 

The corresponding power counting is called the "small scale expansion" (SSE). It is important to 
note that the mass splitting m& — does not vanish in the chiral limit, therefore the explicit in- 
clusion of the delta is a phenomenological extension of CHPT. More recently, it was realized that 
for certain considerations/processes a Lorentz-invariant formulation of baryon chiral perturbation 
theory is advantageous, particular related to the spin sector probed e.g. in doubly virtual Comp- 
ton scattering. A particularly elegant scheme to perform covariant calculations is the so-called 



"infrared regularization" (IR) of [17]. In [57] a consistent extension of the infrared regularization 
method in the presence of spin-3/2 was given. It was in particular shown that the contribution of the 
non-propagating spin- 1/2 components of the Rarita-Schwinger field can be completely absorbed in 
the polynomial terms stemming from the most general effective chiral Lagrangian (see also the 
recent work in [^8|]). 



Here, I report on some results obtained in Q59|]. In that paper, the nucleon and the delta mass 



as well as the corresponding sigma terms were analyzed to fourth order in the small parameter e 
(for related early work on this problem see [|(j]). These explicit representations of and m& 
serve as chiral extrapolation functions to analyze the lattice data from the MILC collaboration p6| ] 
for unphysical pion masses as low as M n ~ 350 MeV (note that we treat the SU(3) data as if they 



were SU(2) - the same comments as in Sec. 3.1 apply). Clearly, such extrapolation functions based 



on chiral perturbation theory cease to make sense at too large values of the quark (pion) masses, 
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but as we will demonstrate later, we can nicely capture the trend of these data. Before proceeding, 
it should be stressed that in the theory with explicit spin-3/2 degrees of freedom one has more 
oeprators and thus more LECs than in nucleon CHPT - similar to the case when going from SU(2) 
to SU(3) in CHPT discussed earlier. 

The corresponding expansions to £?(s 4 ) for the nucleon and the delta mass take the form [59] 



/// 



m 



- Ac x M\ - AmIAqDx - D 2 A 3 - Ae x (X )M 4 - £1 A 4 , - 4M*$E 2 + m N - loop + m A - loop , 



m A = ml - Aa x M\ - AM^AqD\ 



n A A 3 



Ae\(X)Ml ~ E i A o ~ 4M*A$E$ + + m 



N— loop . A— loop 

i r _L_ 111 ' 



(3.14) 



where I have made explicit the contribution from the various contact interactions. Consider first 



the nucleon. Its representation looks very similar to the one given in Eq. (p. 1|). There are, however, 
some important differences. First, the appearance of the delta in the loops also generates additional 
renormalizations ~ Aq (n > 1) (accompanied by combinations of dimension three (Di 2,Df 2 an d 
dimension four (Ei2,Ef 2 ) LECs) that can be absorbed in the LECs of the theory without delta (this 
is discussed in detail in [joT]]) . These LECs can be chosen in such away that the delta contribution 
to the nucleon mass starts at ^(M 4 ). Second, since there is an explicit delta-loop contribution, the 
values of the LECs c, differ from the ones determined earlier by the explicit delta contribution that 



has been worked out in [20]. Third, the delta loops introduce new parameters, in particular the 
leading axial nucleon-delta coupling called ca and the dimension two couplings b\ and from the 
effective nNA Lagrangian appear in the expression for rrifj, utilizing 



= c a{ wt wi n Vn + h.c. } 



5? 



(2) 



TlNA 



'3 iwL f + — iw'uv iL>V + • • 
r mo 



(3.15) 
(3.16) 



with w^j = (t'u^}/2 and w^ v = (x' [D^,m v ])/2, where t', i = 1,2,3 denote the Pauli matrices in 
isospin space. Here, denotes the spin-3/2 field (for more details, see [p0|]). The corresponding 
fourth-order diagrams for the delta self-energy are shown in the left panel of Fig. |5| In particular, 
one has dimension two and four LECs corresponding to the ones in the nucleon mass case, for 
completeness I display the pertinent terms of the dimension two effective Lagrangian, 



_5f 



(2) 



7l& 



ai(X+)- £2 {(u a up)D a DP +h.c.} + |( M 2 ) + . . . 



(3.17) 

so that the a, of the delta correspond to the c, of the nucleon and similarly e\ is a combination of 
dimension four LECs as it is the case for e\ discussed earlier. 

We are now in the position to analyze the pion mass dependence of the nucleon and delta mass 
given in Eqs. ( 3.14 , 3T4| ). They contain a certain number of LECs, some of which are (not very 



accurately) known from the study of pion-nucleon scattering in the heavy baryon SSE Q62[|. In 
addition to the known values at the physical point we take the data from MILC ] [2"6| ] for the nucleon 
and the delta as function of the pion mass and try to describe these with LECs of natural size. Such 
a description is indeed possible, as shown in Fig. |5[ So we do not intend least-square fits here but 
rather try to find out whether the existing data shown in this figure can be consistently described by 
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Figure 5: Left panel: Fourth order contribution to the delta mass. Shown are the contact interactions, pion 
loops (dashed lines) with intermediate nucleons (single lines) and intermediate deltas (double lines). The 
relevant LECs are also shown. Right panel: The nucleon mass (black line) and the (real part of the) delta 
mass (red line) as a function of the pion mass. The filled diamonds denote their physical values at the 
physical pion mass. The dashed line is the chiral extrapolation for the A based on SU(6) as explained in 
the text. The filled squares and circles are the MILC data [p6|. The filled triangles are the recent data from 
QCDSF [§3[]. 



our mass formulas with LECs of natural size. We stress again that a more refined analysis of e.g. 
pion-nucleon scattering in the covariant SSE is mandatory to put stringent constraints on certain 
combinations of the LECs. The parameters corresponding to these curves are: (i) the NA mass 
splitting in the chiral limit, Ao = 0.33 GeV, indicating a slightly larger NA mass splitting in the 
chiral limit than at the physical point, (ii) The LECs from the pion-nucleon Lagrangian are given 
by c\ = — 0.8GeV _1 (which is within the uncertainty of the values determined in e.g. [18]) and 
e\ = C2 = 0, C3 = 0.5 GeV -1 . The small values of C23 are consistent with resonance saturation 
studies of [20] and the fits in [62]. The fourth-order LEC e\ induces the largest uncertainty - even 
a small value of e\ leads to a sizeable contribution at larger pion masses, (iii) The three axial NA 
LECs are ca = 1-1, £3 = 0.75 GeV -1 and b$ = —0.75. (iv) The axial AA coupling is found to be 
\\a = 2, which is not far from the SU(6) or large-A^c value Ha = 9^/5 = 2.28. Furthermore we get 
e\ = -lGeV~ 3 , a x = -0.3 GeV" 1 , A/4 + B = a 2 /4 + a 3 + ... = 0.5 GeV" 1 . These are all natural 
values. It is interesting to note that a\ is markedly smaller than c\, although both couplings should 
be equal in the SU(6) limit. Still, it is interesting to study the strict SU(6) limit. In that case, one 
would have a\ = c\ = — 0.8GeV _1 and Ha = 2.28. As can be seen from the dashed line in Fig. ||, 
the assumption of strict SU(6) symmetry is clearly at odds with the MILC data, indicating that a\ 
and c\ indeed seem to have different values. Also shown in Fig. || are the recent QCDSF data for 
ntA, which were not used in the fit but are nicely consistent with our extrapolation function. Note 
also that the QCDSF data are based on two-flavor simulations and are not very different from the 
MILC data in the region of overlap. This further supports our assumption on the treatment of the 
MILC data. We stress again that the resulting values of the LECs are to be considered indicative and 
a more detailed analysis employing also constraints from other physical processes should follow. 
From the small value of the LEC a\ one immediately deduces that the 71A sigma term appears to be 
significantly smaller than its nucleon cousin because at leading order in the quark mass expansion 
we have g k n = —Ac\M\ + ... and o n ^ = —Aa\M\ + .... It is clear that this interesting observation 
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deserves further study. Finally we note that the sigma term for the nucleon resulting from this 
"rough" fit is found as 

c KN = 48.9 MeV , (3.18) 
to order e A . We note that this value is consistent with the classical result of Ref. [Q|, which was 



confirmed in [18] in a heavy baryon CHPT analysis of pion-nucleon scattering and in [E9f] in a 



CHPT analysis of lattice data. It is also in agreement with the recent CHPT analysis of the three- 



flavor MILC data, see [3J. For the 7iA sigma term we get c n & = 20. 6 MeV. Note that model studies 



give slightly larger values between 28 and 35 Mev [|65|, |66p. Again, these results need to be refined 
and bolstered by more detailed precise fits to the lattice data including also error and correlation 
analysis including also lattice data on other observables - the mass data alone are not sufficient to 
precisely pin down all parameters. Such an analysis, however, is not yet available. Note also that 
chiral extrapolation functions for the baryon octet and decuplet masses in partially quenched QCD 



based on an EFT with deltas are discussed in [p7L 68]. 



4. Quark mass dependence of nuclear forces 

Because of the smallness of the up and down quark masses, one does not expect significant 
changes in systems of pions or pions and one nucleon when the quark masses are set to zero (with 
the exception of well understood chiral singularities like e.g. in the pion radius or the nucleon 
polarizabilities). The situation is more complicated for systems of two (or more) nucleons. Here, 
I report on some work [|^] that is mostly concerned with the properties of the deuteron and the 
S-wave scattering lengths as a function of the quark (pion) mass. These questions are not only of 
academic interest, but also of practical use for interpolating results from lattice gauge theory. E.g. 
the S-wave scattering lengths have been calculated on the lattice using the quenched approximation 



[]70p . Another interesting application is related to imposing bounds on the time-dependence of some 
fundamental coupling constants from the two-nucleon (NN) sector, as discussed in [f7l|]. 
To address these issues, consider the chiral two-nucleon potential at next-to-leading order (NLO). 
It is given in terms of one- and two-pion exchanges (OPE and TPE, respectively) and four-nucleon 
contact interactions 

V = V PE + VtPE + ^contact • (4. 1) 

In this approach, we have to deal with explicit and implicit quark mass dependences. This can be 
most easily understood by looking at the OPE and contact interactions: 

V PE = \—7$2+---dl6 + --- Ti-T 2 77 ( 4 -2) 

Contact = Cs + MlDs + {C T +M 2 n D T ) (di • a 2 ) + . . . (4.3) 

The OPE exhibits both types of quark mass dependences. The pion propagator becomes Coulomb- 
like in the chiral limit, l/(q 2 + M 2 ) — > l/q 2 . The implicit pion mass dependence enters at NLO 
through the pion-nucleon coupling constant (note that the quark mass dependence of the nucleon 
mass only enters at NNLO) expressed through the pion mass dependence of g& /F n in terms of the 
quantity 
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Figure 6: Deuteron BE versus the pion mass. The shaded areas show the allowed values. The light shaded 
band corresponds to our main result with the uncertainty due to the unknown LECs Dj j. The dark shaded 
band gives the additional uncertainty due to the uncertainty of dxf,. The heavy dot shows the BE for the 
physical case M n = M K . 



Here I4 and die are LECs related to pion and pion-nucleon interactions, and the value of the varying 
pion mass is denoted by M n in order to distinguish it from the physical one denoted by M n . In 



particular, d\(, has been determined in various fits to describe kN — > kkN data, see Q54J] . The 
four-nucleon contact terms exhibit an explicit quark mass dependence, These NLO M\ corrections 
to the leading operators ~ Cs,Cj are parameterized by the LECs D$ and Dj, respectively. These 
LECs can at present only be estimated using dimensional analysis and resonance saturation J72|]. 
For the explicit expressions of the TPE, see 



It is now straightforward to solve the Lippmann-Schwinger equation for the bound and scattering 



states utilizing the properly regularized chiral EFT potential, Eq. ( |4.2| ). The deuteron binding 
energy (BE) as a function of the pion mass is shown in Fig|| We find that the deuteron is stronger 
bound in the chiral limit than in the real world, 

= 9.6 ± 1.9^j q MeV , (4.5) 

where the first indicated error refers to the uncertainty in the value of Z)3 S| and d\e being set to 
its average value while the second indicated error shows the additional uncertainty due to the un- 
certainty in the determination of d\(,. Note that in [|73|] a larger variation of the LECs D$j was 
considered, leading to the possibility that the deuteron can become unbound in the chiral limit. 
Also, for massless pions, the higher partial waves are no longer suppressed by the centrifugal bar- 
rier, 8i(k) ~ k (I = 1,2, . . .) due to the Coulomb-like pion propagator. Consequently, one could also 
have bound states in partial waves other than Si — 3 Di (the deuteron channel). However, we did 
not find other bound states. Last but not least, we found smaller (in magnitude) and more natural 
values for the two S-wave scattering lengths in the chiral limit, 

. +0.0 s +0.2 

«c.i.( 1 5o) = -4.1±1.6_ 04 fm, a c ,.( 3 5 1 ) = 1.5±0.4_ 03 fm. (4.6) 
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The pion mass dependence of the scattering lengths and their inverse is shown in Fig. |7[ As can be 
seen from that figure, one needs lattice data for pion masses below 300 MeV to perform a stable 
interpolation to the physical value of M K . We conclude that nuclear physics in the chiral limit is 
much more natural than in the real world. 




M n [MeV] M n [MeV] 

Figure 7: Left panel: The S-wave scattering lengths as functions of pion mass. The shaded areas represent 
the allowed values according to our analysis. The light shaded band corresponds to our main result with 
die = —1.23 GeV 2 and the uncertainty due to the unknown LECs £>s.t- The dark shaded band gives the 
uncertainty if, in addition to variation of As,r, the LEC die is varied in the range from did — —0.91 GeV -2 
to did = —1-76 GeV -2 . The heavy dots corresponds to the values in the real world. The Uiangles refer to 
lattice QCD results from [|70|], Left panel: The S-wave inverse scattering lengths as functions of the pion 
mass. Notation as before. 



5. Summary and outlook 

In this talk, I have discussed the foundations and some applications of the quark mass depen- 
dence of baryon properties, based on chiral perturbation theory and extensions thereof. This can be 
summarized as follows: 

• Baryon CHPT is a mature field in the light quark sector. For the light quark flavors up and 
down, the explicit symmetry breaking is generally small at the physical point. This leads to 
the observation that fourth-order, one-loop calculations provide accurate and unambiguous 
extrapolation functions for lattice QCD. In the three flavor sector, the situation is less satis- 
factory, nonetheless a variety of observables can be studied without resorting to resummation 
techniques, coupled channel analysis or models. The formalism for a covariant inclusion of 
matter fields for spin- 1/2 and spin-3/2 fields exist. Because of decoupling, only the ground 
state baryons can contribute to the leading non-analytic terms in the chiral expansion. If 
one explicitely includes the spin-3/2 decuplet, one must therefore respect decoupling, as 
discussed in Sec. [l[ 
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It is important to realize that observables are linked by general operator structures. Stated 
differently, the low-energy constants (LECs) are universal, they contribute to different pro- 
cesses. For that reason and given the present scarcity of low pion mass lattice simulations, 
one should perform global fits to observables combined with input from phenomenology. 
The quark mass dependence of any observable is given in terms of loop contributions and 
local operators parameterized by class II LECs, the so-called symmetry breakers. The dy- 
namical (class I) LECs corresponding to operators that do not vanish in the chiral limit can 
best be determined by phenomenological analysis. I have tried to stress this issue by consid- 
ering the determination of some of the dimension two pion-nucleon LECs, their appearance 
in the nucleon mass and in loop corrections to the nucleon isovector magnetic moment. When 
considering three flavor CHPT, it is important to realize that the SU(3) LECs are related to 
the SU(2) ones by matching conditions - this allows to get some constraints when analyzing 
baryon octet observables. 

Given the existing lattice data and our present knowledge of the LECs and the quark mass 
dependence of important input parameters like the pion decay constant or the axial-vector 
coupling constant, I conclude that CHPT gives chiral extrapolation functions with a small 
and moderate theoretical error for M n < 400 and 500 MeV, respectively. However, in certain 
cases a smaller window of stability arises at a given order in the chiral expansion, see e.g. 



the discussion of gA in Sec. 3.3. One of the strengths of the effective field theory approach is 
the possibility to systematically work out theoretical uncertainties - almost all figures shown 
here underline this important issue. 

Chiral nuclear effective field theory allows to analyze the quark mass dependence of few- 
nucleon systems, some results for the two-nucleon sector were discussed in Sec. @. The 
two-nucleon force appears to be more natural in the chiral limit than for physical quark 
masses. Also, the observation that both S-wave scattering lengths vanish simultaneously for 
M n ~ 175 MeV has led to the speculation of an QCD IR limit cycle in the three-nucleon 



system [|74|]. Its relevance for lattice QCD was already stressed by Ken Wilson at Lat- 



tice 2004 [75]. 



All this can be summarized in one simple sentence: To connect lattice QCD with real QCD 
(nature) in well-defined and precise manner, we need lattice data at low pion masses - most 
studies performed so far (including the ones presented here) are merely of pioneering nature, 
eventually paving the way for truly ab initio calculations of baryon properties. 
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